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FOREWORD 





A summary report is presented on work performed during the period 
January to September 1974 on Contract NAS8-29877, "Float-Zone Processing 
in a Weightless Environment." The reported work set forth the results 
of evaluations to determine: 1) some practical limits on the size of 

Si crystals that can be grown in a Skylab-type space vehicle, and 
2) the economic impact of large, space-produced Si crystals on the 
electronics Industry. 

This study program is sponsored by the George C. Marshall Space Flight 
Center, National Aeronautics and Space Administration, Huntsville, 
Alabama. Mr. M. Davidson is the COR director of the study. Dr. A. A. 
Fowle, Dr. J. S. Haggerty, Mr. P. F. Strong and Dr. G. Rudenberg of 
Arthur D. Little, Inc., and Prof. R. Kronauer of Harvard University 
are the investigators. 
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1.0 SUMMARY 


1.1 PURPOSE AND SCOPE 

Tlie objective of this work is to evaluate and maximize tiie potential 
benefits of producing faultless silicon crystals in a space environment 
by the float-zone process. A perceived benefit Is the possibility of 
producing larger diameter crystals than can be acliieved on earth whlcli 
may translate into less expensive crystal elements for final use. By 
incorporating relative rotation between the crystal and feed stock, 
forced convective currents may be induced i i the melt which can result 
in a radial distribution of axial heat flux at the freezing Interface 
which reduces grown-in stresses for crystals of the larger sizes that are 
possible to produce in the near zero gravity environment. Specifically 
the objectives of the program were: 

1. Analyze the stability of floating zones to determine the 
practical diameters of Si crystals that can be processed 
in a near zero gravity environment. These diameters will 
be related to rotation rates based on estimates set by 
stability criteria. 

2. For the diameters determined in (1), carry out an economic 
analysis to project the ten year impact to the electronics 
industry. This analysis will include the expected advan- 
tages on the Individual manufacturing phases of integrated 
circuits. 

1.2 CONCLUSIONS AND RECOmENDATIONS 


1.2,1 The Technology of Float-Zone Crystal Growing Process in Space 

Stable molten zones in practical configurations for crystals as large 
as 0.304 m (12 inch) radius can be maintained with limited but useful 
rotation rates in a float-zone process in a near zero gravity environ- 
ment. 

The use of isorotation or counterrotation, rotational rates, and crystal 
and feed stock diameter ratios have been examined qualitatively and seml- 
quantitatively for their influence on producing stress-free crystals of 
uniform composition. Some tentative specification of these process 
variables for best results have been made; however, a combination of 
experimentation and more analysis is needed before a specification of 
these variables for improved and/or larger space grown crystals can be 
made with confidence. 

As a first priority, we recommend some experiments based on model studies 
in order to provide further information relating to the stability of and 
circulation within the molten zone of a silicon crystal process in space. 
The control of fluid motion within the molten zone is critical to a 



successful process and, therefore, a complete imderstanding of tills 
motion Is necessary to any further, useful .-inalyscs directed to the 
optimization of the process. These experiments are recommended as 
the most effective next step In a program leading to tlie design of 
a space experiment. 

1.2.2 Economic Impact 

There is considerable value in large diameter uniform crystals and wafers 
for the electronics Industry. This Is already evident by the shift from 
5 cm (2 inch) to 7.6 cm (3 Inch) diameter crystals by the industry. The 
present study shows that a number of applications can be identified where 
much larger wafers with uniform characteristics have commercial value. 

The principal value of the large wafers is derived from lower chip costs 
because of reduced per unit area processing costs. 

The analysis of the crystal growth process through its various stages 
from polycrystalline feed stock to polished wafers shows that the 
advantage of increased crystal diameters is questionable when based 
strictly on wafer area per unit length of boule or unit weight of poly- 
crystalline feed stock. When made long enough, a float-zone process can 
be considerably more efficient than a Czochralskl process in converting 
polycrystalline feed to acceptably doped single crystals. Offsetting 
this advantage are Increased wafer thicknesses and kerf losses with large 
diameter crystals. The principal savings or value derived from using 
large diameter wafers stems from reduced device processing costs per 
unit area. Secondarily, the efficiency of utilizing area improves with 
larger wafer diameters. 

The economic Impact study was based on substituting 15.2 cm (6 inch) 
wafers for the currently used 7.6 cm(3 inch) wafers. If one-third 
of the total number of smaller wafers were replaced by an equivalent 
area of larger wafers, the annual impact would be about $250-$290 million 
and i\ cumulative ten-year im- ct would be $6.7-$7.3 billion. The use 
of larger, more complex devices will reduce the cost of finished electronic 
packages. The economic Impact from reduced finished package costs is 
extremely tenuous; however, based on past experience it could equal 
the savings derived from reduced chip processing costs. 
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Tliere are several advantages that are potentially realizable for a 
floating-zone crystal growth process carried out under near zero gravity 
conditions. First, there are the factors that arc inherent to a floating 
zone process whether gravity forces are present or not. One of the prln- i 

cipal advantages is that higher purities arc possible witli co-'.tainerloss 
crystal growth processes from high temperature melts. Another is ""he 
Inherently uniform axial dopant levels which result from a zone leveling * 

mode of growth. Second, there are factors arising uniquely as a con- 
sequence of the near weightless environment of a space laboratory. One 
is that buoyancy- induced convection is effectively eliminated from the 
melt. Lack of control over this transport mechanism is responsible for 
many common crystallographic defects. Another is that the growth of 
crystals larger in diameter than can be processed on earth can be enter- * 

tained because of the absence of the destabilizing effects of gravity 
forces. The availability of larger crystals can potentially be trans- 
lated to economic advantage through the value, per se, of large crystal 
elements or by virtue of higher yields of acceptable crystal chips from 
the boule. 

In order to meet the goal of producing large crystals, it is necessary 
to overcome the deleterious effects of thermally- Induced, residual 
stresses in the crystal that become emphasized with scale. An ideal 
solution to this problem is achieved if the conditions of the process 
are arranged to produce a freezing interface that is a plane perpendi- 
cular to the axis of growth. 

In a first generation space laboratory for processing large Si crystals, 
the use of an incandescent heat source is envisioned. This source will 
heat the melt via radiation and, as the silicon melt is opaque to the 
emitted radiation, the heating will take plcCe at the surface. In tn.e 
absence of rotation, heat transport by solid conduction will convey heat 
to interior regions of the melt but this mechanism is inadequate for 
large diameter crystals. 

Rotation of the crystal and feed stock provides a means for effecting j 

"in-depth" heating of the melt zone. Rotation also reduces the azimuthal i 

gradients that would exist in the abse.ice of rotation. Re' ^bive angular ^ 

rotational rates between the feed and crystal and/or different . stal 
and feed stock diameters result in circulating currents within ti». melt 
zone and attendant forced convective heat transport. The technological 
task that presents itself is to define the rotational rates and geometric 
characteristics of the crystal and feed stock that will result in the 
circulation currents that are most favorable for producing a stress-free 
crystal that is also free of crystallographic defects. Although rotation 
may be used to advantage to induce convective heat transport to the 
inner regions of the melt, it also results in destabilizing centrifugal 
forces and can produce periodic variations in the flow field within the 
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melt zone. The type and degree of rotation must, thetefore, be con- 
trolled so that the stability limits of the melt are not exceeded and 
the spatial and temporal variations of velocity associated with the 
forced circulation currents within the melt are attenuated to the 
degree necessary to prevent variations in the solute concentration 
within the growing crystal. 
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3.0 TECHNICAL INVESTIGATION 


3.1 APPROACH 


The approach followed In the execution of the technical invest iK^tl Ion Ik 
outlined below: 

• Carry out dimensional analyses to determine order of 
magnitude Influences of the forces and heat transport 
mechanisms active in tiie melt zone. ConaCruct the Inde- 
pendent dimensionless groups which sets limits on the 
stability of the melt and characterize the velocity and 
temperature fields within a stable melt. 

• Review the literature relating to the stability of 
rotating and non-rotating liquid columns. Extend the 
existing analyses o Include an Important mode of 
instability hitherto unreported. Establish the numeri- 
cal values for the dimensionless groups governing the 
stability of the melt zone. 

• Review the literature relating to the nature of the flow 
In the melt resulting from rotation. Establish the dif- 
ferential equations which determine the velocity fields 
in the melt. As the complexity of a complete solution 
precludes anything but an extensive program involving 
numerical analysis and computers, carry out simple 
limiting/type analyses and use the insights provided by 
the literature dealing with experiments and theory rela- 
tive to similar flows to gain the best possible under- 
standing of the flow within the molten zone anu the 
important physical factors that influence this flow. 

• Establish the differential equations which determine the 
temperature field within the melt. These equations are 
relatively easy to sca’ e if the velocity field is estab- 
lished. In lieu of a means for calculating the real 
velocity field, carry out heat transfer analyses of some 
simple flow models of the melt to gain understanding of 
the relative Importance of convective and conductive heat 
transfer mechanisms, temperature distributions, etc. 

• On the basis of the preceding, postulate the type of 
rotation needed to Induce the most favorable circulation 
currents in the melt. 

• Suggest necessary and promising avenues for Investigation 
in future work 



3 .2 DIMENSIONAL ANALYSIS 


Dimensional analysis has been applied to establish the independent 
dimensionless groups that may govern the behavior of the float-zone 
process in a near weightless environment in order to: (1) provide 

order-of -magnitude insight into the relative magnitude of the 
forces responsible for the absolute stability of the melt zone and 
the circulation within a stable melt; (2) provide order-of-maguitude 
insights into the heat transfer mechanisms responsible for the temper- 
ature field within a stable melt; and (3) determine the similitude 
relationships required for scale model experiments. 

The technical support for the results follow from a straightforward 
application of dimensional analysis common to engineering practice. 

The pertinent results follow. 

3.2.1 Nomenclature 


a > radius: m 


C^(T) * temperature dependent specific heat of the melt 

at constant pressure: j-Kg”^-K" ‘ 

h • height of the melt zone: m 

g - local gravity constant: N-Kg ^ 

k(T) ■ temperature dependent thermal conductivity 
coefficient: watt-m ^-K ^ 


L 


N 


B 





q 


II 


(■ 



■ latent heat of fusion of the melt: J - Kg” 

* Bond number: dimensionless 

” Prandtl number: dimensionless 

" Reynolds number: diuiensionless 

“ Weber number: dimensionless 

>• incident radiation from the source as distributed 
in wavelength and over the melt surface as a 

function of radius ratio, and axial distance 

z -2 ^ 

ratio, e • wAtt - m 
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r 

T 

T 

o 

AT 

V 


z 


reflection coefficient for the melt as dcpeiulent 
on the wavelength X of t^e emitted radiation 
from the source, the local surface temperature 
of the melt T, and the angle of incidence of the 
oncoming radiation expressed as a function of 
2 

7- ; dimension’ ;ss 
h 

radius: m 

temperature : K 

melting or freezing temperature: K 

maximum temperature diff rence between any two 
radial locations within the melt: K 

velocity nnccclated with the circulation within 
the melt: m - sec ^ 

dimension measured in axial direction: m 


a(X.T) 


B 

e(T) 


Y 

X 

P(T) 

0(T) 


absorption coefficient for the melt as dependent 
on the wavelength X of the emitted radiation from 
thr source and the local surface temperature T of 
the melt: dimensionless 

volume coefficient of expansion of the melt: K ^ 


total hemispherical emisslvity coefficient 
(integrated spectral value) as dependent on the 
local surface temperature T of the melt: 
dimensionless 

-2 -4 

Stefan-Boltzmann constant: watt-m -K 

wavelength: m 

temperature dependent density of the melt: Kg~m 

temperature dependent surface tension coefficient 
N-m"^ 


y(T) 


Subscripts : 


temperature dependent coefflci,.it of viscosity: 

„ -1 -1 
Kg-m -sec 


1 refers to the crystal 

2 refers to the feed stock 
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3.2.2 Application to Stability Requirements 

The stability of the Belt tone Is determined solely by the dimensionless 
parameters : 


Weber number, 


P 


This result is subject to four restrictions: (1) gravity forces are 

small In respect to surface tension forces; (2) buoyant forces are 
small In respect to surface tension forces; (3) viscous forces atten- 
dant on the circulation Induced In the melt by rotation are small in 
respect to surface tension forces; and (4) the radial and axial veloci- 
ties associated with the flow within the melt are small in respect to 
the peripheral speed of rotation. These restrictions are met If 


Bond number, N 


Pgtj^ h 


« 1 


pgB&Tr, 


« 1 


yri 


« 1 


'l “l 


« 1 


Assuming that stability Is insured for a Weber number and the geometric 
*^2 h 

ratios — and — of tha order unity, we will show in following sections 
*“l 'l 

that these restrictions are satisfied for SI crystals of large diameter 
processed in a gravity field appropriate to an earth orbiting space 
laboratory. 
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3«2.3 Application to Fluid Circulation Patterns 

An Incandescent heat source applies no applicable forces of electromag- 
netic origin to the molten zone. In an operational regime approaching 
the limits of melt stability, the Influence of buoyant and gravity 
forces have been assessed as unimportant. The effects of variations In 
the surface tension coefficient on fluid motion are also neglected. 
Accordingly, the circulation of the melt Is determined by Inertia and 
viscous forces and, therefore, the circulation patterns within the 
melt are determined by: 


Reynolds number, 


3.2.4 Application to the Temperature Field 

The temperature field within a stable, axlsymmetrlc melt In vacuo which 
Is opaque to the Incident radiation from an Incandescent source Is 
determined by: 


Reynolds number, N_ 


pri 


Prandtl number, 
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Similitude in respect to the stability of a melt requires for both the 
model and prototype that 




3 2 

P 0.^ 


= constant 


(The same value for both model and prototype) 

2 


pgr, 


N, 


B 


P g B AT rj^ 


« 1 


« 1 


^ ^1 “l 
o 

V 


« 1 


ri 0.^^ 


« 1 


^"2 

— “ constant 


^1 


constant 
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constant 


Similitude in respect to the velocity field within a stable molt 
requires for both the model and prototype that 


p r^ 0,^ 


= constant 


constant 


= constant 


= constant 


Similitude in respect to the temperature field within a stable melt 
requires for both the model and prototype that 


a q r 


= constant 


c Y T r 
o o 


= constant 


P V L r, 
1 . 


= constant 


P r, w, 

N„ = = constant 

R y 


constant 


constant 


11 


h 


constant 


* constant 


3.3 STABILITY ANALYSIS 

(D* 

In our previous work for NASA , the absolute stability limits for 
molten zones of various shapes were discussed in detail but the 
destabilizing effects of rotation were not a consideration. The 
results of the previous work are reproduced in Figure 1 which shows 
the family of zone shapes which obtain under steady pulling conditions 
in zero gravity. The geometric limit for absolute stability is indi- 
cated by the dotted circle. The circle represents the stability limit 
for a non-rotating system where surface tension forces dominate, and, 
if exceeded, the molten zone may pinch off into two separate caps. 

( 2 ) 

Hocking discusses the case of a rotating cylindrical column of 
liquid and derives the stability criteria for axisymmetric longitudinal 
disturbances as well as for disturbances that are the same in all 
planes perpendicular to the axis of the cylinder that leave the cen- 
troid of the perturbed shape on the axis. In the latter case he 
obtains the curious theoretical result that a region of neutral 
staMllty exists in the case of an Inviscld liquid which joins the 
• gion of Instability in case the viscosity is finite. For real 
liquids his analysis leads to the same condition that we obtained by 
the simpler means presented in Appendix A, namely that the column is 
stable if 



o 


3 2 

pa (i3 


> 


1 

3 


whei. s 


Tj^ ’* a ^ radius of zone 

= u )2 * “ angular velocity of zone 

In this case the cross section of the molten zone is perturbed to a 
near elliptical shape; that is, by the lowest order radial disturbance. 

Hocking's result for axisymmetric disturbances applies to infinitely 
long columns of liquid. In the case of columns of finite length, the 
disturbance must be a standing wave which does not alter the volume of 
the molten zone so the wavelength must meet the condition 


* 


Superscript numbers in parentheses refer to Section 5.0, REFERENCES 
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5 6 


Maximum Stable 
Zone Diameter 


Limit of Absolute 
Zone Stability 


Maximum Zone 
Height When 
Fiber and Feed 
Rod Have Equal 
Diameters 


Maximum Stable 
Zone Height 


Feed Rod 


FIGURE 1 CALCULATED EQUILIBRIUM ZONE SHAPES PLOTTED IN 
TERMS OF THE ZONE HEIGHT (hi AND RADIUS (r) 
NORMALIZED TO THE FIBER RADIUS. THE CONTOURS 


ZONE HEIGHTS BETWEEN THE FIBER AND THE FEED ROD. 


nX 


h 


where 


n is Integer 

X is the wavelength 

h is the height of the zone 

Then the condition for stability becomes, for n = 1, 


3 2 

pa (ii 



This condition for stability was tested in the experiments of Carruthers 
and Grasso^^^, and they showed it to be applicable over a useful range 
of their experiments. 


In the present program we have examined theoretically the stability of 
the column with respect to another type of axial distrubance, hitherto 
unreported in the literature. In this analysis we consider an unsym- 
metrlc disturbance corresponding to the condition 



= h 


This is a longer wavelength disturbance than any possible axisyimnetrlc 
disturbance and must be asymmetric on account of the requirement to 
preserve the volume of the melt. In this case the molten zone is per- 
turbed to a shape like that of a beam with hinged ends whirling In its 
first bending mode. In Appendix A the condition for stability Is shown 
to be 


3 2 

pa u) 



Graphs of the two stability limits under conditions of rotation are 
shown in Figure 2. It is to be noted that the requirement 


must be met for stability even in the absence of rotation. This is the 
classical result of Rayleigh and is also Indicated by the singularity 
in Rocking's criterion. Moreover, for absolute stability, as opposed 
to possible metastable solutions, our previous work set the limit 

— < 4.3 for a cylindrical molten zone. However, as Figure 2 shows, the 
a 

asymmetric disturbance leads to instability for lower rotation rates 
over the range of h/a of most practical interest. 
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From the above, we might choose a condition of — * n as a practical 

a 

working number to satisfy requirements for the in-dcpth heating of the 
melt and stability under conditions of rotation. In this case, the 
least stable limit of those investigated (L, of Figure 2) reduces to 

a^ 0) 

— ^ * 1, or a Weber number limit of unity. It is to be noted tiuit 

these stability analyses, although providing useful quantitative cri- 
teria for engineering purposes, are limited. They do not account for 
the effects of Internal circulation within the melt or for a variety 
of possible ratios of the rotational rates and the radii of crystal 
and feedstock. To account properly for these effects is a formidable 
analytical problem. It is a problem where the most efficacious solu- 
tion is most likely derived from tests. Some suggested experimental 
methods are described in Section 3.6. 

3.4 STABILITY LIMITED MAXIMUM CRYSTAL ROTATION 


On the basis of the arguments, limitations, and results presented in 
Sections 3.2 and 3.3, a calculation of the maximum rotational rate vs. 
crystal radius has been made for Si molten zones in a radius range of 
0.038 m (1.5 in) to 0.304 m (12 in) with results appearing in Table 1. 
The basis for these calculations are: 


3 2 

_ pa h) 


'W 


1 


N 


B 


0 


r 


1 


a 


— ( 1)2 * to 


h 

a 


It 


o = 7.2 X 10 "^ N/m 


g 


V “ 
P 

9.8 X 10 ^ m/sec^ 


2 X 10 ^ Kg/m - sec 
■ 2.2 X 10^ Kg/m^ 

s 10 ^ times value at earth's surface 
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3 10~^ K~^ 

AT 10^ K 

3.5 HEAT TRANSFER AND FLUID FLOW ANALYSIS 


3.5.1 General Approach 

Heat transfer analysis has two major purposes in application to the 
float-zone processing of large SI crystals. First, It provides the 
means to evaluate various design approaches for producing the In-depth 
heating of the molten zone that is necessary to produce a stress-free 
crystal and, thereby, best methods may be selected. Second, thermal 
analysis must be applied to design a process system that is minimized 
In respect to size, weight and power consumption. In this program 
thermal analysis Is addressed to the first purpose only. Rigorous 
application to the second purpose follows logically sometime after the 
practicality of the space processing of large Si crystals is demonstrated 
to hold real promise. Some estimates of the process system requirements 
have been made by MSFC, and they are regarded to be consistent with the 
capabilities of a next generation Sky-Lab. 

The Ideal conditions of heat transfer to be sought are those that result 
in a freezing Interface that is a plane perpendicular to the axis of 
growth, for it is this condition that minimizes the thermally induced 
stresses in the growing crystal. To achieve this ideal requires that 
the heat flux plus the change of phase heat release per unit area at 
all local areas of the freezing interface match that transferred to the 
growing crystal under a unique circumstance; i.e., that the freezing 
Interface is a plane perpendicular to the axis. It is clear that the 
solution to this problem involves the heat transfer characteristics 
of the whole system consisting of the polycrystalline feed stock, the 
molten zone, and the crystal. 

Solution to the heat transfer problem associated with the feed stock 
and crystal is fairly straightforward. As these materials, at their 
operating temperatures, are opaque to the TR radiation emitted from the 
molten zone, the solution involves phonon conduction in an axisymmetric 
system in vacuo with radiation boundary conditions at the exposed 
cylindrical surfaces* and with either a heat flux or a fixed temperature 
boundary condition at the ends. We, among others, have developed com- 
puter programs to solve this problem routinely. 

The difficult problem is the analysis of the heat transfer within the 
melt. With an opaque melt, we have heat transfer by solid conduction 
and forced convective heat transfer induced by the rotation of the mol- 
ten zone. The importance of the heat transport contribution via con- 
vection requires that the velocity field within the melt be known, and, 
more to the point, be manipulated to achieve the desired condition for 
producing stress-free crystals of uniform composition. 


Other boundary conditions can be easily accommodated. 
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TABLE 1 

PARAMETER LIMITS FOR STABLE, ROTATING, CYLINDRICAL SI MELTS 


a 

m 

u 

rad /sec 

atu 

m/sec 



2 

pgPATa /a 

liaui/o 

0.038 

2.55 

0.097 

3640 

0.0014 

4.32 X 10"® 

2.70 X 10"^ 

0.076 

0.91 

0.069 

5190 

0.0055 

1.73 X 10"^ 

1.92 X 10"^ 

0.152 

0.32 

0.049 

7230 

0.0218 

6.91 X lO"^ 

1.36 X lO"^ 

0.304 

0.12 

0.036 

10700 

0.0870 

2.76 X 10"® 

1.0 X 10"^ 

An 

examination 

of the 

last three 

columns 

of Table 1 shows 

that three of 


the four restrictions necessary to the assumption that inertia and surface 
tension forces determine stability are satisfied. Satisfaction of the 

V 

fourth restriction, i.e., — <<1, is questionable. This is discussed in 

ou) 

Section 3.5.4 following. The maximum speeds of rotation appearing in Table 1 
set the limits within which favorable circulation currents must be sought. 
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If the Influence of temperature on the properties which govern the 
circulation within the melt be ignored (as we believe it can be), then, 
the solution to the problem of determining the velocity field within 
the molten zone can be approached independently of the solution of the 
thermal problem. In this approach, the material properties which 
govern the circulation within the melt are fixed at their values mea- 
sured at the melting temperature. Finally, fluid mechanical considera- 
tions applied to the flow within a rotating, stable melt show the flow 
to take place in the laminar regime which Introduces some measure of 
simplicity In an otherwise very complex problem. 

3.5.2 Analysis of Temperature Field within the Molten Zone 

The differential equations governing the combined convective and con- 
ductive heat transport in a rotating, opaque axlsymmetrlc molten zone 
circulating In laminar flow are developed in Appendix B. The equation 
satisfied by the temperature is 

2 

^ J. ar ^ ■ »' s “1 L ■ ar ^ y ■ jr] . „ 

[a r'^ r' 3r' 3z'^J ^ L 3r’ ^ 3z J 


where 


V/r^ 


subscript 


r refers to radial direction 
z refers to axial direction 

and all other symbols are defined in Section 3.2.1. 


The solution to Equation 1 subject to the appropriate boundary condi- 
tions (see Appendix B) can be gained by application of numerical tech- 
niques and computer calculation provided that the velocity field 

II I * 

(V and V as a function of r and z ) are known. The solution is a 
r z 

If 

description of the temperatures in the melt as a function of r and z , 
or the temperature field. 
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I 


4 


It is of some Interest to note that the dimensionless coefficient 
2 

p C ti u) 

^ is the product of the Reynolds number and the Prandtl 

k 

number, a result impossible to forecast by dimensional analysis. How- 
ever, as the Reynolds number is an independent determinant of the 
velocity field, the Reynolds number and Prandtl nuuber are independent 
determinants of the temperature field as set forth in beet ion 3.2.4. 


Two factors complicate the solution of Equation 1. F 
ous, the contours of the exposed surface of the moltt 
rotating conditions should be known. Second, there 
which specified the isothermal melting and freezing it. 
coordinates and meets the heat flux boundary conditions at these inter- 
faces. It will require iteration in the computational process in order 
to determine this unique solution. 


, to be rigor- 
e under 

nique solution 
. ce In space 


The development above was carried out to reveal what is involved in, 
and methods for, achieving a complete solution to the thermal behavior 
of the molten zone. Actual solution was beyond the scope of this pro- 
gram. 


3.5.3 Evaluation of Azimuthal Temperature Gradients 

As shown in Section 3.4, the stability limited maximum allowable rate 
of crystal and feed-stock rotation varies monotonlcally from 2.55 rad/ 
sec (24.35 rpm) for a radius of 0.038 m (1.5 in) to 0.12 rac/sec 
(1.15 rpm) for a radius of 0.304 m (12 in). As a result of concern 
that such low rotational rates might Impose the need for careful con- 
trol of azimuthal temperature gradients, an analysis of this problem 
was completed. It was found in ^.11 cases that even these low rotational 
rates result In a considerable reduction (smaller by a factor of 20 to 
30) of the temperature gradient that would exist if no rotation were 
present. The analysis and detailed results appear in Appendix C. 

3.5.4 Analysis of Flow Field within the Molten Zone 

The problem here is to define the flow field within the molten zone 
that is most suitable for producing stress-free crystals of uniform 
composition. A review of the pertinent literature* shows that there are 
no analytical or experimental results that are adequate to this defini- 
tion; however, an interpretation of this literature provides useful 
guidance. 

One important simplifying characteristic of the flow field of interest 
is that it is laminar. Bo^h experiment and theory support this conten- 
tion. In fact, the literature shows that the Reynold's numbers asso- 
ciated with the fluid currents within the molten zone are of order Vv 


*See REFERENCES 
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2 

P *^2^ 

with N_ = — as usual. That is, the radial and axial cotnpont'ntK 

VI 

of the velocities Induced In the by rotation times the effective 

hydraulic radius of the current cross sections in loc..l rcRions is 
small in respect to the product of the peripheral speed of tlie crystal, 
r^ times the crystal radius. As shown In Section 3.4, the maximum 

value of Nj^ set by stability for crystal sizes in a range of Interest 

vary from about 3600 to 10,000. Therefore, the Keynolds numbers asso- 
ciated with the radial and axial currents within the molten zone are of 
order 60 to 100. The experience with steady, incompressible chann.l 
flows characterized by Reynolds numbers in this range shows them to be 
laminar. However, just because the currents flow In :he laminar range 
does not mean they are simple in pattern. In fact, they are very com- 
plex. An idea of this complexity is given by the experiments of Hide 
and Titman'^^) , Carruthers and Nassau(^) and Carruthers and Grasso(3), 
even though none of these experiments truly duplicate the conditions 
imposed on the molten zone in a typical float zone process. 

The physical basis for and general character of the flow patterns within 
the molten zone may be described as follows. Differences in the radii 
or angular velocity of rotation of either the crystal or feed stock give 
rise to differences in the radial pressure gradient appearing next to 
the solld-l:*quid Interfaces resulting from the tangential velocities 
Induced near these surfaces by viscous shear. The differences in the 
pressure that results near these boundaries at any radius give.s rise to 
axial , assure gradients, and axial flows take place In response to 
these r .adientb. The general picture of these flows in an axial plane 
through its axis of syimnetry is one of rising and falling axial flows 
separated radially by cylindrical surfaces where viscous shear is great 
and separated axially by radial plane where viscous shear is also large. 
The flow regions of high shear next to the freezing and melting inter- 
faces are commonly called the Ekman boundary layers. The axial flow 
next to the outer surface of the melt is called the Rossby layer by 
Hide and Tltman. In cases where the rotational rates of the crystal and 
feed stock are different, another radial plane of high shear regior is 
located Intermediate between the melting and freezing interfaces. Of 
course, a tangential velocity distribution Is present throughout the 
flow field having local velocities that match those of the crystal and 
feed stock at the melting and freezing interfaces. The picture in an 
axial plane, then, is one of rotating cells separated by axial and 
radial boundary ayers. The cells are ofte.a called Tay lor-Proudman 
cells In recognition of the pioneering work of these scientists on 
viscous flows of the type under discussion. 

The flow picture described Is qualitatively correct but Plmpllfied. 
Depending on the Reynolds number characteristics of thj flow, tne 
boundary layers an become large (at low values of Reynolds number) in 
respect to the linear dimensions of the molten zone and the cellular 
and boundary layer divisions become smeared. Moreover, at high 
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Reynolds numbers, shear waves can develop in regions of high viscous 
shear, as observed by Hide and Titman, which fact further complicates 
the flow picture. 

With the qualitative picture of the behavior of the flow field in mind, 
the question becomes, how can one select the rotational and dimensional 
characteristics of the system to best achieve in-depth heating in order 
to minimize thermal stresses in the growing crystals while obviating 
variations in solute concentration Imparted by flow fluctuations? In 
simple terms the objective is to convey the heat absorbed at the exter- 
nal surfaces of the melt to interior regions while suppressing temporal 
and spatial variations in the velocity field. 

Figure 3 illustrates In schematic form some more promising arrangements 
for achieving the desired result. They appear in the order of their 
expected merits. The expected qualitative features of the fluid flow and 
heat transfer characteristics of the system shown in Figure 3a is 
presented below as an example of what might be made to take place in a 
desirable design. 

Starting at a location near the outer surface of the molten zone next to 
the freezing interface, we note an axial flow directed toward an inter- 
mediate radial plane of symmetry where the tangential component of 
vel'^city is zero everywhere. This axial stream is heated by the absorbed 
radiation from the source, and, by virtue of convective heat transport, 
a major portion of the net radiation absorbed is transferred to the 
intermediate ]>lane of symmetry. Here the axial stream diverts in direc- 
tion to flow radially inward in a boundary layer that can be likened to 
a typical Ekman layer. The radial flow in this boundary layer peels 
off as it proceeds toward the axis of symmetry giving rise to a uniform 
axial flow proceeding in the direction of the freezing Interface. As 
this radial flow is bounded by an adiabatic surface at z = 0, the tem- 
perature of this flow is unchanged. Within the axial flow proceeding 
from the intermediate plane of symmetry to the freezing interface, we 
note that the ratio of heat transport by forced convection to solid 
conduction reduces, becoming zero at the interface. The net effect is 
that the enthalpy flux entering radially at the Intermediate plane, which 
is nearly the same as the absorbed radiation, is uniformly distributed 
over the freezing interface where it is released by heat transfer via 
solid conduction. To complete the picture, the axial flow approaching 
the freezing Interface is diverted to flow radially outward in the 
Ekman boundary layer next to this interface, thereafter to turn and flow 
again axially next to the exterior surface of the melt. 

In this idealized nicture, the intermediate radial plane of symmetry 
is isothermal wit . temperature elevated above the melting point. The 
total heat absorbed in the outer layer of the melt is distributed uni- 
formly over the freezing interface. If the crystal is insulated to pre- 
vent radial gradients, then, the freezing interface will be a plane 
perpendicular to the axis of symmetry — the desired result. With the 
reasonable assumption that the heat flux contribution due to the latent 
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FIGURE 3 SCHEMATIC REPRESENTATION OF AXIAL AND RADIAL 
COMPONENTS OF VELOCITY FIELD IN MOLTEN ZONE 
FOR SOME MORE PROMISING FLOAT-ZONE SYSTEMS 
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heat of solidification is relatively small, the temperature and flow 
fields in the molten zone between the intermediate plane and the 
melting interface is symmetrical with respect to those described. 

The qualitative nature of the velocity fields associated with the other 
systems depicted in Figures 3a and 3b can be constructed with reference 
to that given for Figure 3a. In Figure 3b the path length for solid 
conduction from the outer surface of the melt to the freezing interface 
is reduced to take possible advantage of this fact for in-depth heating. 
Figure 3c is presented ns a hedge against the possibility that the 
location of the intermediate planes that divide the cells in Figures 3a 
and 3b may fluctuate in space about some mean location, with consequent 
variations in solute appearing in the growing crystal. The idea 
advanced here is that the system of Figure 3c is probably more stable 
in this regard although its heat transport characteristics are less 
desirable. 

The description of the flow field given above borrows heavily from the 
results of the prior experimenters cited and on the supporting analysis 
of Hide and Titman. However, it is to be noted that the boundary condi- 
tions on the exterior arfaces of the equivalent molten zone in these 
prior experiments deviate significantly from the real cases of interest 
in the present investigation. 

As a final observation, it is noted that the radial and axial components 
of velocity within the molten zone as Indicated by the measurements of 
Carruthers et al are small in respect to the peripheral velocity of the 

crystal. That is, ^ << 1, a state satisfying the conditions for the 

stability criterion advanced in Section 3.4. However, against this, the 
work of Hide and Titman predict that the radial and axial components of 
velocity in local regions within the melt can approach the peripheral 
speed. This being the case, local distortion of the outer surface will 
result from these local dynamic effects, and the stability criterion 
might be modified somewhat. Whether these effects on stability are 
significant is a debatable point. We believe that they are probably 
not significant, but this contention needs to be tested by experiments 
of the type proposed in Section 3.6. 

An analytical description of the velocity field within the molten zone 
is provided by the appropriate form of the Navier-Stokes equations 
together with the equation of continuity. These equations are developed 
in Appendix D. 

For an incompressible, steady, laminar, axisyrametric flow with material 
properties independent of temperature, the Navier-Stokes equations com- 
bined with the continuity equation in dimensionless form are 
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The equation of continuity alone is 

A U V , . 

3r ' / 3z 




( 2 ) 


We note here the appearance of the Reynold's number, 


pri 0.^ 


, a result 


anticipated by the dimensional analysis presented in Section 3.2.3. 


Few analytical solutions to Equations la-lc and 2 are available. Those 
that are involve relatively very simple flow systems where the 
constraints reduce the active variables to manageable proportions. In 
principle, numerical techniques involving machine computation will serve 
to achieve a solution for our cases of interest. The introduction of 
the concepts of a stream function and potential function as an aid to 
solution does not apply in our case for the potential function does not 
exist in a system having vorticity, although the stream function exists. 
Numerical solution to problems of the type we are concerned with utilize 
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relaxation teclmiques in botli time ami space coordinates. TIh* problem 
statement treats the system in unsteady state (we have to include terms 
3V 

of the form — ) and relaxes it to tliu steady-state solution employing 
o C 

tlie appropriate initial and boundary conditions. Tlie paper by Ors/ag 
and Israell^^^ reviews tlie work in tills area. 

Upon review of the literatcic pertaininR to solution of similar fluid 
flow problems, one is impressed with the complex nature of possible 
flows which may result. The work of Pearsonv?) on the viscous tlow 
between two rotating coaxial disks having inf.'nile radii makes this 
point. Even for this relatively simple system which has relatively 
straightforward characteristic solutions to the Navier-Stokes equations, 
Pearson shows equilibrium solutions are highly depende t on the 
Reynolds number based on the disk spacing which is the only geometric 
dimension in the system. Moreover, the equilibrium solution in some 
cases, characterizied by a higher Reynolds number, was shown to depend 
on the start-up condition; that is, which disk is assumed to Influence 
the flow first. This analytical problem is idealized, but the result 
does give an indication of the sensitivity of the flow pattern to slight 
variations in the speeds of rotation of the disk. 

A determination of the velocity field within the molten zone by numeri- 
cal solution of the governing equations is outside the scope of the 
current program. The work presented was carried out with the objective 
of gaining maximum insight into the nature of flow field short of that 
provided by computer solution of the governing equations and of defining 
the magnitude of the problem to be faced should computer solution be 
attempted . 

3.6 FUTURE WORK 


In order to pursue the technical investigation beyond what has been 
accomplished in this program requires: (1) a more accurate determina- 

tion of the stability of the molten zone under the variety of configu- 
rations and rotational rates that may have utility; and (2) a more 
accurate determination of the velocity field within the molten zones 
associated with each promising system. Emphasis should be placed on 
achieving the latter objective. 

The work accomplished clearly shows the magnitude of the difficulty 
associated with a "quantum- jump" in the sophistication of analysis 
represented by computer solution to the Navier-Stokes equations. 
Confronted with this difficulty, it would be easy to attempt solutions 
to problems without adequate knowledge whether the results were war- 
ranted or valid. Accordingly, we have concluded th.'t experiments are t 
vital necessity to the definition of useful flow fields in the Si melt. 


26 




» \ 
I 


f 


i 


I 


I 


(3) 

Carruthers and Grasso have carried out useful experiments of the 
type needed. However, these experiments were limited in respect to tlie 
stability conditions that they could examine, and the flow patterns 
that they observed resulted from a boundary condition on the model of 
the molten zone that was significantly different from that which exist 
in a real float-zone process. 

First, as regards stability, these experiments examined the influence 
of perturbations which lead to distorted cross sections with a centroid 
remaining on the axis only. Therefore, they did not examine the 
effects of the type of perturbation that we have determined to be most 
limiting as regards the stability of typical melts of interest. Second, 
these experiments did not examine the Influence of Internal circulation 
on stability. 

An examination of the Carruthers and Grasso experimental results pertain- 
ing to the observed circulation patterns make clear that the experimental 
method did not produce truly the circulation patterns of a real float 
zone. This can be seen by reference to the results presented for equal 
isorotation of equal radii, simulated crystal and feed stock. These 
results showed significant internal circulations of a magnitude roughly 
the same as the case of equal counterrotation. From elementary conside- 
rations, we know that, in the case of a melt with a free surface, no 
circulations of the type described can persist at steady state with 
equal Isoratlon. Anything other than solid body rotation of the complete 
melt would give rise to shear stresses at the interfaces, a net torque 
and resulting rate of change in angular momentum of the melt. The 
Induce currents observed by Carruthers and Grasso are the result of shear 
stresses at the exterior boundary of the simulated melt induced by the 
mineral oil surroundings. It seems reasonable to assume that the 
presence of the oil Influences the circulation patterns in all the cases 
they described. To what degree the patterns they observed are different 
from those of a real float-zone is moot, but their results for equal 
isorotation make the translation in all cases questionable. 

The purpose of experiments we propose in an on-going program would be 
five fold: (1) they would test the results of the stability analyses 

that we have made; (2) they would establish stability limits for system 
configurations and operating conditions not anzlyzed; (3) they would 
test the validity of concepts we have developed as to the nature of the 
circulating currents induced in the melt by various forms of rotation; 

(A) they would provide guidance for the numerical analysis and computer 
solution of the velocity fields within these melts; and (5) they would 
obviate some major limitations of prior experiments reported in the 
literature. 

A brief description of the type of experiments we propose is given in 

Appendix E. The experiments involve the use of scale models and optical ' 

tracing of the motion of observable partfcles suspended in the model 

melt. From our examination of possibilities, we have found that a model 

-3 i 

melt system using water with a simulated crystal radius of 2.7 x 10 m > 

(0.103 in) has useful potential of meeting the objectives stated above. I 
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4.0 ANALYSIS OF ECON(»aC IMPACT 



4.1 SUMMARY 

4.1.1 Purpose and Scope 

We have considered the possible ten-year economic impact of growing larf>c 
float-zone silicon crystals in a gravity-free environment. We describe 
the expected advantages of fabricating large crystals, wafers and chips, 
and then present examples for the most significant crystal diameters, 
characteristics and devices. 

For this analysis, we have only considered the economic impact on devices 
that already have a sizable, developed market, and we have assumed re- 
placement of 5 cm (2 inch) float-zone and 7.6 cm (3 inch) Czochralskl 
crystals by 15.2 cm (6 inch) float-zone crystals. We have not conducted 
an actual volume analysis of specific device types, because the data- 
gatheri- g and reduction would involve an extensive effort. 

The following discussion considers wafer size to be the major influence 
on the ultimate chip cost. A more extensive analysis will be required 
to account for the Influence of space shuttle costs, power costs in 
orbit, the practical (rather than theoretical) yield of good crystals, 
and the actual crystal quality obtainable in a weightless environment. 

The. conclusions presented in our analysis assume that all of the tech- 
nical and economic problems of growing and processing 15.2 cm (6 inch) 
crystals in orbit will be resolved satisfactorily. 

4.1.2 Introduction 


It appears that crystals several times the present 5-cm (2 inch) and 
7.6-cm (3 inch) diameter products currently available could be grown 
in a gravity-free environment. Large diameter, high quality float- 
zone crystals would provide a number of economic advantages for manu- 
facturing semiconductor devices. For example; 

• The float-zone process produces a greater volume of 
acceptable dopant level for any given crystal size. 

• The float-zone process generally produces higher 
quality crystals. 

• Large float-zone crystals could use large diameter 
feedstock, which costs less per kilogram than the 
small diameter feedstock currently used in the float- 
zone process. 
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• The large wafers available from large liiameter crystals 
would have much higher efficiencies of surface uti I i/.aLion, 
which would lead to much greater production of chips per 
unit area - particularly for large scale integrated circuits. 

• The greater nuinber of chips per wafer will greatly reduce the 
cost per chip, because wafer processing costs are essentially 
independent of wafer size. 

• Large, high quality wafers will enable the production of new, 
larger classes of rectifiers, thyristors, and LSI's, and these 
larger, more complex circuits will reduce the overall cost of 
both finished devices and electronics equipment. 

The possibility of growing large float-zone crystals in orbit presents 
at least one major disadvantage that should be considered: 

• The technology necessary to grow large, float-zone crystals 
in an orbital environment is currently beyond the state of 
the art, and considerable research and development effort 
will be required to solve the technical and economic uncer- 
taintier inherent in such processing. 

4.1.3 Method of Economic Analysis 

a. Polycrystalline Feed to Trimmed Boule . As the first step in analyzing 
the value of large float-zone crystals, we have examined the cost of 
converting polysilicon feedstock into both Czochralski and float-zone 
crystals. 

We first consider the cost of producing larger diameter feedstock and 
the possible cost advantage of using such material to produce large 
float-zone crystals versus smaller Czochralski crystals. 


We then consider the cost of producing crystals of equal mass and equal 
diameter and the value per unit mass of feedstock of the trimmed boules. 

Finally, we examine the conversion efficiencies of the two processes, 
with equal size feedstock, per unit mass of trimmed boule. At this 
stage, the output of crystal with acceptable dopant level and resis- 
tivity determines the value of each trimmed boule. 

b. Trlmned Boule to Finished Wafers . The analysis of this manufactur- 
ing stage considers wafer output primarily in terms of surface area per 
unit weigiht of trlmned boule. 

c. Wafer Processing to Finished Chips . This analysis of processing 
discusses the geometry of area utilization, including the output of 
chips per unit area as a function of wafer diameter and the influence 
of ’arger wafers on area utilization for large chips. 
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d. Chip Processing Costs . Ihis manufacturing stage Is discussed In 
terms of device complexity (as a function of the number of processing 
and handling steps ) , device type and wafer size. 

e. Market Impact . Finally, the Influence of producing large float- 
zone crystals In orbit is considered in terms of current and pro'iected 
market size, device cost as a function of wafer size, and the economic 
Impact on selected device types. 

4.1.4 Conclusions 


There is no clear-cut cost advantage for either growing process, because 
there are only marginal savings possible through the different feedstock 
utilization or conversion efficiencies. In fact, the savings inherent 
in each process nearly cancel out by the time each type of crystal has 
been converted to a trimmed boule. 

However, it appears that the number of finished chips available from 
each 15.2 cm (6 inch) wafer could be four tc eight times the number 
of finished chips produced from the currently available 7.6 cm (3 inch) 
wafers. Therefore, because processing costs are very high and are 
essentially independent of wafer size, our analysis leads us to conclude 
that the processing efficiencies made possible by the availability of 

15.2 cm (6 inch) wafers would enable the manufacturers of silicon semi- 
conductor devices to produce such devices for perhaps 12% to 25% of the 
present cost per chip. The theoretical economic impact on the semi- 
conductor Industry of such savings can be extrapolated to a cumulative 
ten-year value of $8 billion to $15 billion. 

4.2 CZOCHRALSKI VERSUS FLOAT- ZONE FOR MASSIVE CRYSTALS 
4.2.1 Introduction 


This section considers the cost of producing crystals from polysilicon 
feedstock, and examines: 

• Feedstock cost as a function of diameter, 

• Crystal cost as a function of feedstock size, and 

• Conversion efficiency to trimmed boule. 

4.2.2 Feedstock Cost 

The cost per unit weight , ' polysilicon feedstock decreases with increasing 
diameter. This results because radial deposition rates are essentially 
Independent of size and the cost is largely determined by capital and 
energy costs for the time spent in the trlchloro-sllane reactor. With 
these deposition kinetics, the weight of silicon produced per unit time 
increases proportionally with diameter. Table 2 summarizes the average 


TABI.E 2 


FEEDSTOCK GROWING TIME AND WEIGHT 
VERSUS FINAL ROD DIAMETER 


Polyslllcon 
Rod Diameter 

Length of 
Deposition Run 

Welglit of 
Silicon Per 
0.3 m (12") Rod 

Average 
Production 
Rate Per 
0.3 m (12") Rod 

cm (inches) 

(days) 

(kg) 

(kg/day) 

5.7 (2.25) 

5.5 

1.8 

0.33 

7.6 (3) 

7 

3.2 

0.46 

10.2 (4) 

12 

5.8 

0.48 

12.7 (5) 

17 

8.5 

0.50 

Source: Arthur D. 

Little, Inc. 
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SI production rates as a function of diameter. The average production 
rate for 12.7-cm (5 Inch) polyslllcon rod Is 50% greater than the rate 
for 5.7-cm (2.25 Inch) rod. 

The 7.6-cm (3 Inch) Czochralskl process currently uses the largest feed 
rod diameter (12.7 cm, 5 Inches), which costs $0.70 to $1.00 per gram, 
while the 5.7 cm (2.25 inch) float-zone process uses the smallest rod, 
which cost $1.30 to $1.50 per gram. Thus, for the earth-based process, 
the float-zone feedstock Is 50% more expensive than the Czochralskl 
feedstock, and the latter crystal-growing process enjoys a cost advantage. 

4.2.3 Crystal Cost 

If both the processes are used to grow 7.6-cm (3 inch) diameter crystals 
of equal length and equal mass, we find that the Czochralskl crystal costs 
less because of lower feedstock cost. However, we must compare the cost 
of growing a 7.6-cm (3 inch Czochralskl crystal against the cost of grow- 
ing a 15.2-cm (6 incii) float-zone crystal in orbit — and only consider the 
growing process Itself. Then, we find that both processes use approxi- 
mately the same diameter feedstock, which would eliminate the feedstock 
cost advantage of the Czochralskl process. 

4.2.4 Conversion Efficiency 

If 7.6-cm (3 Inch) Czochralskl crystals are compared with 15.2-cm (6 Inch) 
float-zone crystals, we must consider \diat percentage of polyslllcon Is 
converted by each growth process to the form of a trimmed Ingot (or boule) 
that has been cropped and ground to a specified diameter and length with 
the desired resistivity for further processing. 

Figures 4 and 5 Illustrate the geometry of crystal and ingot for the 
Czochralskl and float-zone processes, respectively. 

The Czochralskl process is used to produce uniform P-type material and 
both broad-range (2.5:1) and narrow-range (1.7:1) N-type material. When 
all of the crystal-growing and machining losses and recycling of material 
have been accounted for, the weight of acceptable P-type trimmed ingot 
material is about H5% of the initial charge. For the widely used N-type 
crystals, the material with a broad resistivity range provides about 65% 
acceptable trimmed Ingot, while the narrow-range material provides about 
45% of the Initial polyslllcon charge in the form of usable trimmed ingot. 
(The charge weight was assumed to be 7.5 kilograms.) 

Industry experience with 5.7 cm (2.25 inch) float-zone crystals indicates 
that a trimmed ingot provides 75-80% acceptable material from a similar 
initial charge weight (7.5 kg). This compares favorably to the yields 
of the Czochralskl process, depending on the type of crystal grown. 


32 








FIGURE S GEOMETRY OF FLOAT-ZONE CRYSTAL. INGOT AND WAFER 






Because the length of float-zone mate: lal lost to ingot cropping may equal 
two crystal diameters, sm-'ll charge weights would give unacceptable yields 
with 15,2-cm(6 inch) diameter crystals. Figure h sliows tlie t lieon t i ca > 
yield of trimmed ingot available from 7.6-cm (3 inch) .ind 15.2-cm (6 inch) 
float-zone crystals, with charge weiglits ranging from 7 to 100 kilograms, 
tliarges of 50 to 100 kilograms would be required tc increase the ylehl 
from 15.2-cm (6 inch) crystals tt> 75-87% accoptal)le material. 

The 15.2-cm (6 inch) float-zone crystal would provide nearly four times 
the weight of acceptable trimmed ingot per unit length that is available 
from a 7.6-cm (3 inch) crystal — Czochrairki or float-zone. 

4.2.5 Conclusions 


At this point, it would appear that 15.2-cm (6 inch) float-zone crystals 
grown in orbit would be perhaps four times as valuable as 7.6-cm (3 inch) 
Czochralski crystals grown on earth. However, the uncertainties involved 
with growing the 1 rger crystals in orbit may ultimately negate the theo- 
retical advantage that we have just described. Thus, to this stage in 
the manufacturing process, neither process has a clear-cut cost advantage. 
In fact, the average raw material cost is such a small percentage (7-^Z) 
of the cost of a finished device that any savings to this point may have 
only marginal impact on the ultimate cost of finished chips. 

4.3 EFFECTS OF LARGE CRYSTALS ON WAFER FINISH ING 

Industry experience has shown that the number of wafers per unit length 
that can be sliced from a trimmed boule decreases with increasing bouic 
diameter. This results in part from the face that larger diameter wafers 
must be thicker to survive mechanical handling during finishing and chip 
processing. Also, saw thicknesses must be Increased with larger diameters, 
which Increases kerf losses. 

With present technology, approximately eight polished wafers are obtained 
per centimeter (20 per inch; of trimmed ingot length for 7.6-cm (3 inch) 
diameter crystals. We have estimated that approximately five wafers per 
centimeter (12 per inch) could be expected from 15.2-cm (6 i^ich) crystals. 
Wafer thickness was increased from .05 to .08 cm (.018 to .030 inch) and 
saw kerf thickness from .04 to .08 cm (.015 to .030 inch). 

Wafer breakage and material losses occur during sawing, lapping, polishing 
and other wafer finishing procedures. The material losses per unit area 
should not be influenced significantly by wafer size due to the thicker 
wafers to be used. We believe that the increased wafer thickness should 
permit equivalent yields through sawing and lapping (82%) and polishing 
(85%) steps. 
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Percent Yield of Trimmed Ingot 



36 



Tlius, the cost differential between the 7.6-cm (3 inch) diameter wafers 
and the 15.2 cm (6 inch) diameter wafers is reflected in the wafer area 
per unit length and area per laiit mass of trimmed boule. The larger 
wafers generate 2.4 times the area per unit length and 0.6 times the area 
per unit mass of boule. Without specific information about the comparative 
costs of growing 7.6-cm (3 inch) Czochralski crystals on earth and 15.2-em 
(6 inch) float-zone crystals in orbit, it is impossible to make a quanti- 
tative estimate of the relative costs per unit area. It is obvious that 
the larger space-grown crystals will cost more per unit length — probably 
more than the area differential per unit length. Thus, a cost advantage 
or value of large wafers will probably be foiaid In reduced device fabri- 
cation costs on the larger wafers. 

4.4 AREA UTILIZATION BY FINISHED DEVICES 


A geometric analysis of the useful area available from both 7.6-cm (3 inch) 
wafers and 1,5.2-cm (6 inch) wafers — summarized in Table 3 - indicates 

that the larger wafer can enclose from five to eight times the number of 
individual chips than the smaller wafer can. Thus, the value of the larger 
wafer is even greater than the simple 4:1 ratio of surface areas would 
suggest. 

Two factors make the area converted to chips smaller than the nominal area 
of the wafer. A rim, usually C,32 cm (0.125 inch) wide, cannot be used 
because it is damaged during wafer handling. Ihe percentage of rim area 
to total area decreases with increasing wafer diameter. Secondly, square 
or rectangular chips do not fill a circle completely. If any portion of 
a chip intercepts the outer perimeter, it is unusrble. The efficiency 
of utilizing the enclosed area also improves with increasing wafer diameter. 
Th is latter factor is particularly Important for the larger chip sizes, 
which fit inefficiently onto smaller wafers. 

4.5 CHIP PROCESSING COSTS 


Handling and manufacturing costs for wafer processing have proven to be 
almost Independent of wafer size. Thus, not only does the efficiency of 
surface utilization increase with wafer diameter, the processing costs per 
unit area decrease dramatically with larger wafer diameters. 

The cost of processing a wafer or producing a finished de ice is directly 
related to the device complexity. Figure 7 Illustrates .he average 
production cost (including overhead) of wafer processing as a function of 
complexity, expressed by the number of masks needed to process each wafer. 

The linear relationship exists because the major element in the cost per 
’•afer is the labor needed for masking and diffusion transfers. For example, 
the process for producing Metal-Oxide-Semlconductor (MOS) Large-Scale Inte- 
gration (LSI) requires nearly 100 transfers of wafers to accomplish the 
5 to 7 masking steps. Complex bipolar Integrated circuits may require up 
to twice as much labor cost for handling, transfers, loading of processing 
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F:GURE 7 COST OF WAFER PROCESSING VERSUS COMPLEXITY OF PROCESS 
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equipment, inspection and test, because of the 9 to 12 m. sking steps and 
other procedures. The labor concent ranges from 2.5 Co 5 manhours for 
' rocessing to the point of being cut into individual chips. 

As indicated in Figure 7, bare wafers have an average value of approxi- 
mately $6, while processing costs raise the wafer value to $80 to $100. 

It is evident that raw material cost is a very small percentage of total 
Tinished device cost, except for the most simple devices. The p''incipal 
savings or value derived from using large diameter wafers stem from reduced 
processing costs per loiit area, because wafer processing costs have his- 
torrcally proven to be essentially independent of size. It is assumed 
here that wafer processing costs will not Increase markedly in going to 
15.2 cm (6 inch) diameter Secondarily, the efficiency of utilizing area 
Increases with wafer diameter for the reasons cited in Section 4.4. Thus, 
the number of chips produced per wafer increases faster than a simple 
area ratio would predict. Thn cost of a finished chip produced on the 
15.2-cm (6 inch) diameter wafers should be about 1/5 to 1/8 the cost of 
a chip produced on 7.6-cm (3 inch) diameter wafers. 

4.6 MARKET SIZE. GROWTH AND ECONOMIC iMPACT 

The 1974 demand for silicon wafers — valued at about $1.00 per square inch 
(6.5 square cm) — is currently about 200 million square inches (1.3 billion 
square cm). This demand is growing at an average rate of 20Z per year, 
and therefore is expected to Increase sixfold in ten years. Thus, the 
projected cumulative ten-year market is 25 times the current demand or 
about 5 billion square inche«^ (32.3 billici square cm) — a market of about 
$5 billion for raw wafers. 

The 1974 production consists of some 40 mixllon wafers and about 60Z (or 
24 million wafers) were processed to finished devices. At this time, a 
finished processed wafer repre..ents an average saleable value of approxi- 
mately $40, so that processing could represent a total annual value of 
about $1 billion. (A volume analysis for all devices would be necessary 
to establish the actual value addad by processing each year.) 

To estimate the theoretical economic impact of substituting 15.2-cm (6 inch) 
wafers for 7.6-cm (3 inch) wafers, we must assume that the per wafer pro- 
ces.' Ing costs remain essentially the same, and that the larger wafers 
will provide four to eight times the number of finished devices. (Of 
this also assumes equipment capable of processing the larger wafers, and 
the hypothesis of equal processing yields.) The annual economic Impact 
of the 15.2-cm (6 inch) wafers could then be 75-88r« of $1 billion, or 
$750-$880 million. In fact, it is only reasonable to consider replace- 
ment of about one-third of the smaller wafers with the large diameter 
wafers. This would lead to an annual Impact of about $250-$290 million, 
and a cutaulatlve ten-year Impact of $6,7-$7.3 billion. We believe that the 
economic Impact would be realized primarily through reductions in the 
price of various finished devices, rather than in the development of a 
larger market. 
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An effect, for tihlch there is currently no dollar value, is the possibility 
of producing much larger, more complex devices on the larger wafers, which 
will greatly reduce the cost per functica in completed systems and thus 
reduce the overall c st of finished electronic packages. Historically, 
such impact has shown up as 10-40Z reduction in materials and labor costs 
which are generally 25-35Z of finished package costs. At this point, an 
estimate wouIj be rather tenuous, but it could amount to 2-3Z (10Zx25Z). 
Considering that production of semiconductor devices supplies an estimated 
$35 billion in electronic equipment each year, this could result in a 
savings of $700 million to over $1 billion a year. Long-term results over 
the decade could ultimately be $15-25 billion. 

In actuality, such a theoretical savings can be reached only after con- 
siderable research and development to srtisfy all technical and economic 
problems. Currently, all of the crystal growth and wafer finishing pro- 
cedures assumed in this analysis are beyond the state of the art. 
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APPENDIX A 


THE STABILITY OP A ROTATING, CYLINDRICAL MOLTEN ZONE 


In this appendix we consider the stability of cylindrical molten zone 
rotating as a solid body. The zone Is subjected to two forces: surface 

tension, which tends to maintain the cylindrical shape, and centrifugal 
force, which tends to destabilize the zone. We assume that gravity is 
zero or, at least, negligibly small. The method adopted is to consider 
a tiny displacement of the equilibrium surface shape and to see whether 
or not the distortion grcwrs without limit. 


First, we consider the distortion to be a function of radium alone. In 
Figure A1 we show a circle which represents the undistorted cross section 
of the cylinder with the equation (in polar coordinates) 


r = a 


( 1 ) 


along with a cross section of the distorted cylinder represented by 
the equation 


r ■ a + e cos 20 


( 2 ) 


TTie pressure on the undisturbed surface which arises from centrifugal 
force is 


P ” 


2 


U) p 



r dr 


2 2 
po) a 


2 


( 3 ) 


The destabilizing pressure is the change in p due to the change in the 
sh:ipe of the surface and for small changes is given by 


Ap 


2 . 

pu rAr 


or 


Ap * pu a e cos 20 


(A) 

( 5 ) 


Next we must compute the stabilizing torces produced the change of 
curvature of the surface. It is evident from Figure Al that the greatest 
curvature of the disturbed surface is at the greatest radius, so the 
change of curvature does, at least, tend to restore the cross section to 
circular shape. The change in the other principal curvature also produces 
a restoring force dependent on the axial hel^t of zone, but we will neglect 
that in this calculation. 
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The handiest vay to calculate the change In curvature is to use a formula 


1 

R 


+ 2R*^ - rr" 
(r^ ^ r'2)3/2 


( 6 ) 


that gives the curvature of a curve expressed in polar coordinates, 
r > r(6). Primes denote differentiation with respect to o . We have 


r ■ a + c cos 2 0 
r* = -2e sin 26 
r" * ~4e cos 2 9 


and, to first order, 

r^ + 2r'^ - rr” - (l + ~ cos 20 ) 

, 2 ^ '2,3/2 3 3e _ A 

(r + r ) *a (1+"^ co® 2 9 j 


Equations 8 and 9 may be substituted in Equation 6 to give 


1 

R 


( 


1 + M cos 29 ] 


( 


1 + — cos 29 
a 




---f 

R a \ 


1 + — C 03 20) 


Whence 

&(|) - cos 20 . 
a 

At the limit of stability 

oA(-|) - Ap 


(7) 


( 8 ) 

( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


46 



So, from Equations (5) and (12) , 


3 o g cos 2 e 


> pu a c cos 20 


Thus, to ensure stability, the dimensionless group known as the Weber 
number must satisfy the inequality: 


< 3 


Next, we investigate the slttiation which occurs when the cylindrical 
zone of liquid suffers an unsymme rical uisplacement in a direction 
normal to its axis, say in the direction of the x-axls. We take the 
equation of the disturbed cylinder to be: 


2 2 2 
(x - e cos kz) + y ■ a , 


where " h 


and h is the height of the zone.- Figure A2 shows the cylinder with and 
w.'.thout the perturbation as well as the system of axes employed. 

To first oruer Equation (16) is 


2 0 2 2 
X - 2xe cos kz + y ■a 


2 2 
r ~ 2ae cos kz cos 6 ■ a 


whence we find 


2rAr > 2ac cos kz cos 9 


Eqiution (19) may be siibstltuted into Equation (4) to give the destabiliz- 
ing pressure 


4p •• pu at cos kz cos 9 
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As before, the stabilizing pressure is obtained after calculating the 
change in curvature. In this case, we consider the cuivature in a plane 
passing through the z-axis 




213/2 


From Equation (18) , we find 


■ - k e sin kz cos 0 


1,2 . Q 

— j ■ - k e cos kz cos 0 
8z 


1 1-2 . 

— ■ k Ecos kz cos 9 
K 


Since this curvature is zero in the case of the unperturbed cylinder, 
the right hand size of Equation (24) is the change In curvature which 
gives rise to a stabilizing pressure due surface tension. We have, for 
the limit of stability 

2 2 

0 k e cos kz cos 0 ■ pu a c cos kz cos 0 v25 

and hence we obtain the stability criterion for the height of the 
rotating cylindrical melt in terms of the Weber nun^er as follows: 

2 u2 0 

pu) ah ^ _2 


2 3 
po) a 


. 2-a,2 

< . (jj) 
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APPENDIX B 


EQUATIONS GOVERNING THE TEMPERATURE FIELD WITHIN THE MOLTEN ZONE 


The purpose of this appendix is to set forth the equations governins the 
combinsd convective and conductive heat transport in a melt circulating 
in laminar flow. Appropriately cylindrical coordinates arc chosen (r, 

6, z) with: 




velocity components 
thermal conductivity 

rate )f heat generation per unit volume 


temperature 


t = time 


p ■ density 
y “ viscosity 

The basic differential equation governing the thermal behavior of the 
molten zone results from application of the energy equation to an ele- 
mental volume of dimensions dr, rdd, dz. For material properties assumed 
Independent of temperature, the equation satisfied by the temperature is 

f S) ^ ^ ^ ^ fc V « ^ ^ «» « ^ ^ V, T,] 


, II I 3T 

+ q pCp • • • 


For the axlsymmetrlc case, with an opaque melt at steady state. Equation (1) 
reduces to 




0 r ) 


or alternately. 
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r.)^ T 1 3T T 1 PC£ f- ^ ^ „ 8T ''r ^ y 8I + yZi|. 

“'2”* f ^"■72"J ' “ir L’^ ''t 5 ?* r 

L >r 3z -* 


Introducing the condition for continuity of matter 


i (r vp ^ (r V^) - 0 


reduces Equation (2b) to 


I ill + i + i-2. l_££ v -^i-v ^ 


By substituting the non-dimensional variables 


•^1 "l 


where 


*^1 “l 


r^^ ■ radius of crystal 
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( 


■ melting or freezing temperature 
* angular velocity of rotation of crystal 


Equation (4) becomes 


[g2 3^ I 


pCpr^ ***1 
k 


V • ^ + V • ^ 

L r 3r' z 3z' j 


0 (5) 


The boundary conditions for the molten zone are 


" 4 

a q - E o T 
n 


-(If 


sin 0 + — cos 9 

3^ ) 


for the exposed surface of the melt, and 


, 3T I 

- I 2 COS 9 = 


cos 0 + L u p 


11 I 

jr '1 


- I 

3r '2 


T- = T„ = T = T 
12 o 


for the freezlug Interface, and 


k I 2 cos 0 = k 


H ll cos 9 + L u p 


3T I 

Tf '1 


il 1 

3r '2 


T, - T- - T - T 
12 o 
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where 




= the oncoming radiant heat flux per unit area normal to 
the surface of the molten zone. 

0 = angle between the outward directed normal to the bounding 

surface and the positive z direction. 

u =» rate of crystal growth cr pulling rate. 

= radius of feed stock 

( 1)2 = angular velocity of rotation of feed stock- 


subscript 1 refers to the melt side of the interface 

2 refers to the crystal or polycrystal side of the 
interface 

The solution to Equation (5) subject . the boundary conditions (6) through 
(12) can be gained by application of numerical techniques and computer 
calculaticn provided that the velocity field (V^ and Vj as a function of 
r' and z') are known. The solution is a description of the teraperatu’'es 
in the melt as a function of r' and z', or the temperature field. 
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APPENDIX C 


REDUCTION OF AZIMUTHAL TEMPERATURE GRADIENTS BY ROTATION 


We consider the silicon crystal, the molten zone and the feed rod to be 
a long cylinder which is rotating as a solid body. We assume that the 
cylinder is subject to surface heating which is not uniform in azimuth 
and ask what effect rotation has in reducing the azlmuthel temperature 
gradient thus produced. 

The equation satisfied by the temperature T in steady state is 


7 • kVT - C pv 7T * 0 

P 


where 


k is thermal conductivity, watts cm ^ K ^ 

C is heat capacity, joules g~^ 

P -3 

p is density, g csi 

V is velocity vector, cm sec 

The velocity is in the e-dlrection only and has the component 

Vg . u,r 


where 


u) is the angular velocity, radians sec 
r is the radius, cm 


Next, we assume the properties of the material to be constant and write 
Equation 1 in cylindrical coordinates: 


k 3 / 3 t\ ^ k 

r 3r \ 3rj ^ ^2 


3^T - 3T 


It seems simpletit to adopt an Indirect procedure at this juncture in the 
study of the azimuthal gradient problem. We assume a suitable form for 
the temperature distribution and derive the heat input from it by 
computing the radial gradient at the surface. We let 


y(r) e 
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and substitute this expression in Equation 3 to fiid the function y(r) 
which gives the dependence of temperature on the radius. After the 
substi*^ution is made and some obvious algebraic manipulations are per- 
formed, Equation 3 reduces to a form of Bessel's equation: 




1 C pui 


The solution is seen to be 


y = - (■v/n’Kr), 


where 


and is an arbitrary constant (degrees C) which relates to the ampli- 
tude of the gradient. 

The Bessel function cf complex argument appearing in Equation 6 is 
customarily expressed in terms of Kelvin functions. We find the temper- 
ature distribution to be 

T =* T^(berj^ Kr + 1 beij^ Kr) e^® (8) 

■2 

The radial heat flux q" (watts cm ) is given by the formula 


q" = 


3r a 


where a is the radius of the cylinder, cm. 

From Equation 8 we see that 

q" - - ^ \ I ^ 


+ — (ber, Ka + 1 bei, 
a 1 ] 
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I It Is useful CO consider the magnitude of the ratio of the surface 

, temperature « T , to tie heat flux, from Equations 10 and 8, it may be 

seen that 



For short we let the radical appearing in Equation 11 be called u(Ka), 
which function is tabulated in Table Cl. For small values of Ka, the 
function is seen to be equal ^o its argument. For large values of Ka, 
one could appeal to the known properties of the Kelvin functions to 
evaluate u(Ka), hut it is easier to return to Equation 5 and obtain an 
approximation directly. 

First we make the transformation 


1 

2 

2 “ r y 

in Equation 5 to obtain the equation for z: 



( 12 ) 


(13) 


2 

For large values of Kr we may drop the term 3z/4r in Equation 13 to 
obtain the simpler equation 


M - 

dr*^ 


(14) 


whose solution is readily seen to be 


2 • A e 


1±1 

/2 


Kr 


(15) 


where A is an arbitrary con-tant. Therefore, the tetnperature T is given 
by the formula 


i±i 

T - JL, /2 


Kr + 16 


(16) 
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TABLE Cl 


Function to be Used in Determining Azimuthal Temperature Gradients 


Ka 

u(Ka) 

0.0 

0.00000 

0.5 

0.49871 

1.0 

0.96090 

1.5 

1.26132 

2.0 

1.33012 

2.5 

1.27218 

3.0 

1.19586 

3.5 

1.13960 

A.O 

1.10576 

4.5 

1.08689 

5.0 

1.07595 




q 


II 


KAk 


14-1 

vt 


Ka+10 




(17) 


The magnitude of the ratio of surface temperature to heat flux is thus 
found to be 




to first order. 


(18) 


While the results we have derived so far are useful to predict the azi- 
muthal temperature gradient to be expected when the azimuthal distribu- 
tion of heat flux is known or can be estimated, it may be of more value 
in the preliminary stages of syste;.. design to consider the ratio of the 
temperature gradient which appears in the absence of rotation to the 
gradient obtained with rotation ii; the presence of the same azimuthal 
distribution of heat flux. This ratio, which we call the smoothing 
factor f, is given by the relation 


f 


K a 
u(Ka) 


(19) 


For large values of Ka the smoothing factor is given by the approximate 
formula 


f = Ka ^ (20) 

For values of Ka > 4, this approximation produces less than 1% error. 

In Table C2 are shown the values of the sioothing factor f for Ka _< 5. 

In order to evaluate Ka, we use the following values to describe the 
physical properties of silicon: 

k = 0.22 cm"^ k“^ 

p = 2 . 2 g cm ^ 

0 =“ 720 X 10 ^ N cm ^ 

C = 1.^44 J g'^ k"^ 

P 
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and we use the Weber nurabei:, given by the ratio 

3 2 

N = (21) 

w o 

to determine the rotation rate w. Combining Equations 21 and 7, we may 
write 


Ka 


1 

fC \2 

. (o p N, a) 



( 22 ) 


Finally, we take = 1 and find from Equation 22 that 


Ka = 13.7A29 a (a in cm) 


(23) 


For a > 1 cm, it is clear that the smoothing factor can be taken roughly 
to be equal to Ka. Some values of the smoothing factor f and corres- 
ponding rotation rates as a function of crystal radius a are shown in 
Table C3. In spite of the low rotation speeds, the smoothing factors 
are all quite large, a fact which may be attributed to the rather low 
thermal conductivity of the silicon near its melting temperature. 
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TABLE C2 


S oothlng Factor 


a 

f 

0.0 

1.0000 

0.5 

1.0026 

1.0 

1.0407 

1.5 

1.1892 

2.0 

1.5036 

2.5 

1.9651 

3.0 

2.5087 

3.5 

3.0713 

4.0 

3.6174 

4.5 

4.1403 

5.0 

4.6471 


TABLE C3 

Smoothing Factors and Rotational Rates for a Weber Number of Unity 


a, cm 

f 

rpm 

5 

20 

15.5 

10 

24 

5.5 

15 

27 

3.0 

20 

29 

1.9 

25 

30 

1.38 

30 

32 

1.05 
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APPENDIX D 


EQUATIONS GOVERNING THE FLOW FIELD WITHIN THE MOl.TEN ZO NE 


In this appendix we set forth the equations governing the circulation of 
the melt in laminar flow. They are the Navier-Stokes equations and the 
equation of continuity appearing widely in the literature. Appropriately, 
a cylindrical coordinate system (r, 0, z) is chosen with 

Vr* Vq, “ velocity components 

p = density 
U = viscosity 


Axial symmetry, incompressibility am-, steady state is dictated. Gravity 
forces and variations in material properties with temperature are 
assumed to be negligible. Under these conditions, the Navier-Stokes 
equation combined with the continuity equation, is 




DV y" 3V 

_x._§_ + v — " 

3r r ^ 3z 




) 


3r 


+ P 


3^V , 3V V 3^V 

r ^ 1^ r _r i 


P IV 


V 

6 ^ r 0 



^ 1 


’^0 

-f 

3r r 

z 3z/ ~ 

M 

2 t 

3r^ 

m 

3r 

2 

r 


!!Is1 

,z^\ 


(la) 


(lb) 


P 


3V 


3r 


+ V 


av. 

I 

3z 


) 


3z ^ “ 


3^V , 3V 3^V 

^+1 ^ + 


3r 


3r 


The continuity equation alone is 


3z 


(Ic) 




( 2 ) 


Kith the substitution of 


r 

r. 


z 

r. 
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Pti 

We note here the exlstenrc of the Reynolds number, , a result 

u 

anticipated by the dimensional analysis presented In Section 3.2. 

Few analytical solutions to Equations 2 and 3a-3c are available. Those 
that are li volve relatively very simple flow systems where the 
constraints reduce the active variables to manageable proportions. In 
principle, numerical techniques Involving computer computation will 
serve to achieve a solution for our cases of Interest. The Introduc- 
tion of the concents of a stream function and potential function as an 
aid to solution ^ocs not apply In our case for the potential function 
does not exist in a system having vortlcity, although the stream func- 
tion exists. V imerical solution to problems of the type we are concerned 
with utilize relaxation techniques In both time .< i space coordinates. 

The problem statement treats the system in unst .dy state (we have to 

dV 

include terms of the form ) and relaxes it to the steady-ctate solu- 
tion employing the appropriate initia’ and oundary conditions. In our 
case, the initial conditions would be ~ ^^0 " ^^2 ~ ^ everyw.iere except 
.(t the boundaries. 


The bounaary conditions are: 


At the cylindrical surface, r ’ f(z) 


P 0 


(rj + 4) 


where : 

P = static pressure 

0 = the coefficient of purface tension 

and R 2 = the principal radii of curvature of the surface 
of the melt zone under conditions of Insured 
stability. 

At the boundary z ■» 0 (assumed plane) 

V = V = 0 

r z 


At the boundary z = h (assumed plane) 


V = V =0 
r 7 . 


63 


APPENDIX E 


SOME PROPOSED EXPERIMENTS 


In order to model the behavior of a stable Si melt in the spacecraft 
environment, we should Ideally produce an earth-based test system in 
which the Inertial forces due rotation are about equal to the surface 
teriSion forces and each is large in respect to gravity forces. This 
m^ans that the dimensions of the model melt must be relatively very 
small. Backing off from the ideal, we might establish tentatively the 
condition that the inertia, surface tension and gravity forces are all 
the same magnitude in the model melt. This would establish a reasonable 
upper bound on the size and the magnitude of the circulating currents 
that wG can expect to produce in the laboratory. 


For the reasons stated above, 
the conditions: 


N, 


W 


our model melt should meet. 


3 2 

= Weber No. = ^ — — = 1 


as a minimum. 


( 1 ) 


Ng = Bond No. = P S = 1 (2) 

H =1 

a 

2 

N = Reynolds No, = ^ — — as large as possible (3) 

R n 


To meet the condition expressed by the limiting value of the Bond number, 
we wish to select a fluid substance that has the smallest possible ratio 
of p /a in order to experiment, with a model having the largest possible 
dimensions. The Weber nuni>er limit for stability coupled wi h the desire 
for a model melt having the largest possible dimensions also dictates the 
selection of a fluid having a miulraum value of p /o . The desire for a 
large Reynolds number, on the otbof hand, calls for a large value of 
p and a small value of p. On expressing the Reynolds number in terms 
of the Weber nunber, we get 


and with Ny * 1 


‘V " a P O / w 

Nj^ « ^ a p c / p 
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As a result of investigation into substaiices that might best meet the 
conditions above and are, in addition, readily available, easy to use, 
and best be transparent for flow visualization studies, we find that 
water appears as a good candidate. For water at room conditions: 

P = 1 X 10^ Kg/m^ = 1 gm/cm^ 
o = 7,2 X 10 ^ N/m = 72 dyne/cm 

_3 

p = 1 cm. = 1 X 10 Kg/m-sec 

a = \! K ^ = J (1) — = 2.71 X 10~^ ra 
V 2 P8 Y 

= ^a p o /p = A41 

NrP 

u = — 2 ~ 60 rad/sec = 573 rpm 

pa 

A Reynolds number of AOO is down by a factor of about ten from those 
values that can be achieved in stable Si molten zones in a space-based 
process. Examination has revealed that a model system with mercury makes 
possible experiments at a maximum Reynolds number of 2300. Stability 
and gravity limits restrict the radius of the mercury model to some value 
less than about 0.19 cm, and certainly flow visualization studies with 
mercury will be less rewarding. 

However, from what we know now, it is not clear at what value of Reynolds 
number the optimum benefits from rotation is to be achieved. One view 
contends that the optimum lies with larger Reynolds numbers. This view 
is based on the assumption that the ^aracter of the flow would have a 
certain similarity independent of Reynolds number in the Reynold? number 
range of interest. The assumption of the more circulation the better 
may still be true, but as the entire character of the flow may well change 
in the Reynolds number range of interest, this simple assumption is in 
question. In fact, flows characterized by larger Reynolds numbers may 
have shear waves and other undesirable periodic features. It seems more 
reasonable, in the light of our current knowledge, to consider utilizing 
better def.ied and understood flows even though they may not represent 
the optimum that might eventually be achieved. In this light, experiments 
wi^h water still appear to be very useful. 
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